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The hot gas in the IGM produces anisotropics in the Cosmic Microwave Background (CMB) 
through the thermal Sunyaev-Zel'dovich (SZ) effect. The SZ effect is a powerful probe of large-scale 
structure in the universe, and must be carefully subtracted from measurements of the primary CMB 
anisotropics. We use moving-mesh hydrodynamical simulations to study the 3-dimensional statistics 
of the gas, and compute the mean comptonization parameter y and the angular power spectrum 
of the SZ fluctuations, for different cosmologies. We compare these results with predictions using 
the Press-Schechter formalism. We find that the two methods agree approximately, but differ in 
details. We discuss this discrepancy, and show that resolution limits the reliability of our results to 
the 200 < / < 2000 range. For cluster normalized CDM models, we find a mean y-parameter of the 
order of 10 -6 , one order of magnitude below the current observational limits from the COBE/FIRAS 
instrument. For these models, the SZ power spectrum is comparable to the primordial power 
spectrum around I — 2000. It is well below the projected noise for the upcoming MAP satellite, 
and should thus not be a limitation for this mission. It should be easily detectable with the future 
Planck Surveyor mission. We show that groups and filaments (kT < 5 keV) contribute about 50% 
of the SZ power spectrum at I — 500. About half of the SZ power spectrum on these scales are 
produced at redshifts z < 0.1, and can thus be detected and removed using existing catalogs of 
galaxies and X-ray clusters. 



I. INTRODUCTION 



The hot gas in the IGM induces distortions in the 
spectrum of the Cosmic Microwave Background (CMB) 
through inverse compton scattering. This effect, known 
as the thermal Sunyaev-Zel'dovich (SZ) effect |l]J^], is a 
source of secondary anisotropics in the temperature of 
the CMB (see Refs. §-§ for reviews). Because the SZ 
effect is proportional to the integrated pressure of the 
gas, it is a direct probe of the large scale structure in 
the low rcdshift universe. Moreover, it must be carefully 
subtracted from the primary CMB anisotropics, to allow 
the high-precision determination of cosmological param- 
eters with the new generation of CMB experiments (see 
HQ and references therein). 

Thanks to impressive recent observational progress, 
the SZ effect from clusters of galaxies is now well estab- 
lished gg,§-|l§. The statistics of SZ clusters were calcu- 
lated by a number of authors using the Press-Schechter 



(PS) formalism (eg. |ll|-|l5|). Recently, Atrio-Barandera 
& Miicket |l(| used this formalism, along with assump- 
tions about cluster profiles to compute the angular power 
spectrum of the SZ anisotropies for the Einstein-de Sit- 
ter universe. A similar calculation was carried out by 
Komatsu & Kitayama (KK99, hereafter) p7| , who also 
studied the effect of the spatial correlation of clusters and 
cosmological models. 

The statistics of SZ anisotropies have also been studied 
using hydrodynamical simulations. Scaramella et al. [fl8|| , 
and more recently da Silva et al. [|l9| , have used this ap- 
proach to construct SZ maps and study their statistical 
properties. Persi et al. |20) instead used a semi-analytical 
method, consisting of computing the SZ angular power 
spectrum by projecting the 3-dimensional power spec- 
trum of the gas pressure on the sky. 

In this paper, we follow the approach of Persi et al. 
using Moving Mesh Hydrodynamical (MMH) simulations 
pl|p2[ . We focus on the angular power spectrum of the 
SZ effect and study its dependence on cosmology. We 
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compare our results to the Press-Schechter predictions 
derived using the methods of KK99 . We study the red- 
shift dependence of the SZ power spectrum, and estimate 
the contribution of groups and filaments. We also study 
the effect of the finite resolution and finite box size of the 
simulations. Results from projected maps of the SZ effect 
using the same simulations are presented in Seljak et al. 
|p3| . We study the implications of our results for future 
and upcoming CMB missions (see also Refs. p4| , p5[ ). 

This paper is organized as follows. In §||, we briefly 
describe the SZ effect and derive expressions for the inte- 
grated comptonization parameter and the SZ power spec- 
trum. In 



[II 



we describe our different methods used to 
compute these quantities: hydrodynamical simulations, 
the PS formalism, and a simple model with constant bias. 
We present our results in § |TV| , and discuss the limitations 
imposed by the finite resolution and box size of the sim- 
ulations. Our conclusions are summarized in § fv| . 



II. SUNYAEV-ZEL'DOVICH EFFECT 

The SZ effect is produced from the inverse Compton 
scattering of CMB photons (U-§. The resulting change 
in the (thermodynamic) CMB temperature is 



AT 

-TfT = yj(x) 
Jo 



(1) 



where To is the unperturbed CMB temperature, y is the 
comptonization parameter, and j(x) is a spectral func- 
tion defined in terms of x = his/k^To, h is the Planck 
constant and ks is the Boltzmann constant. In the 
nonrelativistic regime, the spectral function is given by 
j(x) = x(e x + l)(e x — 4, which is negative (positive) 
for observation frequencies v below (above) i>a ~ 217 
GHz, for T ~ 2.725 K. In the Rayleigh- Jeans (RJ) limit 
(x <C 1), j{x) ~ —2. The comptonization parameter is 
given by 



V 



dl 



m e c 2 



m e c 2 



dl p e 



(2) 



where <jt is the Thomson cross-section, n e , T e and p e 
are the number density, temperature and thermal pres- 
sure of the electrons, respectively, and the integral is over 
the physical line-of-sight distance dl. 

We consider a general FRW background cosmology 
with a scale parameter defined as a = R/Rq, where 
R is the scale radius at time t and Rq is its present 
value. The Friedmann equation implies that da — 

1/2 

Ho (l — + Slma -1 + fUa 2 ) dt where O = U m + Qji, 
f2 m , and are the present total, matter, and vacuum 
density in units of the critical density p c = 3H 2 /(8irG). 
As usual, the Hubble constant today is parametrized by 
H = 100 h km s _1 Mpc -1 . It is related to the present 
scale radius by Ro = c/(kHq), where n 2 = 1 — 0, 1, and 
il — 1 in a open, flat, and closed cosmology, respectively. 



The comoving distance x, the conformal time r, the light 
travel time t, and the physical distance I are then related 
by dl = cdt = cadr = ad\- 

With these conventions, and assuming that the elec- 
trons and ions are in thermal equilibrium, equation (j^) 
becomes 



y = a T \ adx 



P 



~ p e m p Tn e c 



(3) 



where p is the gas mass density, T is the gas temperature, 
and p~ l = n e J (p/m p )is the number of electrons per pro- 
ton mass. Equation (|3|) can be written in the convenient 
form 



(4) 



y = yo dxT p a 



where T p = pT j~p is the gas density- weighted tempera- 
ture, and p = Pc^bO, 3 - The overbar denotes a spatial 
average and Qb is the present baryon density parameter. 
The constant y Q is given by 



2/o 



p, e m p m e c 2 



1.710 x 10" 



-16 



0.05 



1.136 

Pe 



KT 1 Mpc" 1 , 



(5) 



where the central value for p e was chosen to correspond 
to a He fraction by mass of 0.24, and that for £1/, to agree 
with Big Bang Nucleosynthesis constraints 

The mean comptonization parameter y can be directly 
measured from the distortion of the CMB spectrum (see 
Ref. p7l for a review), and is given by 



y = Vo dxT p a 



(6) 



It can thus be computed directly from the history of 
the volume-averaged density- weighted temperature T p . 
The gas in groups and filaments is at a temperature of 
the order of 10 7 K (or 1 keV), and thus induce a y- 
parameter of the order of 10~ 6 over a cosmological dis- 
tance of cHq 1 ~ 3000 h' 1 Mpc (see Eq. [§). This is 
one order of magnitude below the current upper limit of 
y < 1.5 X 10" 5 (95% CL) from the COBE/FIRAS instru- 
ment |||. 

The CMB temperature fluctuations produced by the 
SZ effect are quantified by their spherical harmon- 
ics coefficients a; m , which are defined by AT(n) = 
T -1 ^2 lm ai m Yi m (n). The angular power spectrum of the 
SZ effect is then C; = (|a; TO | 2 ), where the brackets de- 
note an ensemble average. Since most of the SZ fluc- 
tuations occur on small angular scales, we can use the 
small angle approximation and consider the Fourier co- 
efficients AT (1) = / rf 2 nAT(n)e' l n They are related to 

the power spectrum by (AT(1) AT*(1')) ^ T 2 {2tt) 2 5 {2 ^ (1- 
\')Ci, where 8^ denotes the 2-dimensional Dirac-delta 
function. The SZ temperature variance is then a 2 -. = 



2 



{(AT/To) J = Ei(2i + l)Ci/(47r) ~ fdllCi/(2w). Since, 

as we will see, T p a~ 2 varies slowly in cosmic time scale 
and since the pressure fluctuations occur on scales much 
smaller than the horizon scale, we can apply Limber's 
equation in Fourier space (eg. Ref. to equation (^) 
and obtain, 



C l ^j\x)yt / d X T p P p - >X a- 4 r 



(7) 



where r = i?o sinh(x-R ( 7 1 ), X> an d -Ro sin(xi?^" 1 ) are the 
comoving angular diameter distances in an open, flat, 
and closed cosmology, respectively, and P p (k, X ) is the 3- 
dimensional power spectrum of the pressure fluctuations, 
at a given comoving distance x- I n general, we define the 
3-dimensional power spectrum P q {k) of a quantity q by 



<«r g (k)C(k')> = (2nfd^(k-k')P q (k), 



(8) 



where S q (k) = J d 3 X(5 ? (x)e 4k x , and S q = (q — q)/q. 
With these conventions, the variance is cj^ = (6%) = 
J d 3 kP q (k)/(2ir) 3 . For aflat universe, equation (0) agrees 
with the expression of Persi et al. p3]. The SZ power 
spectrum can thus be readily computed from the history 
of the mean density- weighted temperature T p (x) an d of 
the pressure power spectrum P p (k,x)- 



per grid cell. The code that optimized for parallel pro- 
cessing, which is straightforward due to the regular mesh 
structure. The moving mesh provides linear compression 
factors of about 10, which correspond to compression fac- 
tors of about 10 3 in density. Note that this code does not 
include the effects of cooling and feedback of the gas. 

We ran three simulation with 128 3 curvilinear cells, 
corresponding to erg -normalized SCDM, ACDM, and 
OCDM models. The simulation parameters are listed in 
table |. Note that in all cases, the shape parameter for 
the linear power spectrum was set to T — Vt m h p6|p5[ . 
The simulation output was saved at z — 0,0.5,1,2,4,8 
and 16, and was used to compute 3-dimensional statis- 
tics. 

To test the resolution of the simulation, we compared 
the power spectrum of the dark matter density fluctua- 
tions P p DM(k) (defined in Eq. || with q = pdm) from 
the simulations to that from the Peacock & Dodds |3(| 
fitting formula. The results for the ACDM are shown on 
figure ^, and are similar for the other three models. The 
simulation power spectrum agrees well with the fitting 
formula for 0.2 < k < 2 h Mpc -1 at all redshifts. For 
k < 0.2 and k > 2 h Mpc -1 , the simulations are limited 
by the finite size of the box and the finite resolution, re- 
spectively. We will use these limits below, to study the 
effect of these limitations on the SZ power spectrum. 



III. METHODS 



B. Press— Schechter Formalism 



A. Simulations 

We used the MMH code written by Pen [flHH], which 
was developed by merging concepts from earlier hydro- 
dynamic methods. Grid-based algorithms feature low 
computational cost and high resolution per grid element, 
but have difficulties providing the large dynamic range 
in length scales necessary for cosmological applications. 
On the other hand, particle-based schemes, such as the 
Smooth Particle Hydrodynamics (for a review see Ref. 
pfj) fix their resolution in mass elements rather than 
in space and are able to resolve dense regions. How- 
ever, due to the development of shear and vorticity, the 
nearest neighbors of particles change in time and must 
be determined dynamically at each time step at a large 
computational cost. 

To resolve these problems, several approaches have re- 
cently been suggested |3^-j3§] . The MMH code combines 
the advantages of both the particle and grid-based ap- 
proaches by deforming a grid mesh along potential flow 
lines. It provides a twenty fold increase in resolution over 
previous Cartesian grid Eulerian schemes, while main- 
taining regular grid conditions everywhere |22| . The grid 
is structured in a way that allows the use of high resolu- 
tion shock capturing TVD schemes (see for example Ref. 
|p3 and references therein.) at a low computational cost 



It is useful to compare the simulation results with an- 
alytic calculations based on the Press-Schechter (PS) 
formalism |57|. We compute the angular power spec- 
trum and the mean Comptonization parameter, using the 
methods of KK99 and Barbosa et al. Jl^], respectively. 
For definitiveness, we adopt the spherical isothermal j3 
model with the gaussian-like filter for the gas density 
distribution in a cluster, 



Pgas(0 = PgasO 



-3/3/2 



(9) 



where R and r c are the virial radius and the core radius 
of a cluster, respectively, and a fudge factor £ = 4/n is 
taken to properly normalize the gas mass enclosed in a 
cluster jTi] |. We employed a self-similar model for the 
cluster evolution (3£|. Note that other evolution mod- 
els yield spectra that differ only at small angular scales 
(I > 2000) JIt]]. The gas mass fraction of objects is taken 
to be the cosmological mean, i.e., Q(,/f2 m . 

The volume-averaged density-weighted temperature is 
given by 



1 



M„ 



dMM MM > z) T(M,z), (10) 
dM K ' 
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where ~p = 2.775 f2 OT h 2 M Q Mpc -3 is the present mean 
mass density of the universe, dn/dM is the PS mass func- 
tion which gives the comoving number density of col- 
lapsed objects of mass M at z. T is computed by the 
virial temperature given by 



k B T(M,z) = 5.2 /T 1 
x(l 



A c (z) 



1/3 



18tt 2 
z) a 1 / 3 keV, 



M 



10 15 h- 1 Mr. 



2/3 



(11) 



where A c (z) is the mean mass density of a collapsed ob- 
ject at z in units of p £! m (l + z) 3 While Barbosa 
et al. [|l5| used (3 ~ 5/6, we adopt (3 = 2/3 according to 
KK99. 

The limits M m [ n and M max should be taken to fit the 
resolved mass range in the simulation. The mass enclosed 
in the spherical top-hat filter with comoving wavenumbcr 
k is 



M 



An 



-fH) 



\k 



3.6 x 10 13 



1 h Mpc" 



©• 



(12) 



Since the fc-range of confidence in the simulation is ap- 
proximately 0.2 < k < 2 h Mpc" 1 (see §[IIA), equation 
(H) gives M min ~ 4.5 x 10 12 Q m h~ x M Q and M max ~ 
4.5 x 10 15 il m hr x Mq. This mass range is used for cal- 
culating the angular power spectrum, the mean Comp- 
tonization parameter, and the density-weighted tempera- 
ture. A more detailed inspection of Figure || reveals that 
the resolution of the simulations depends on redshift, and 
involves a power law cutoff in k rather than a sharp cut- 
off. This must be ke pt in mind when comparing the two 
methods (see §[VD). 



C. Constant Bias Model 



It is also useful to consider a simple model with con- 
stant bias. The bias b p of the pressure with respect to 
the DM density can be defined as 



b 2 Jk,z) 



^ P p (k,z) 

PpDM (^! Z) 



(13) 



and generally depends both on wave number k and red- 
shift z. In this simple model, we assume that b p is inde- 
pendent of both k and z, and replace the pressure power 
spectrum P p (k,z) in Equation ((?]), by b p P p unj(k, z), 
where P p dm is evaluated using the Peacock & Dodds 
fitting formula |pq| . This has the advantage of allowing 
us to extend the contribution to the SZ power s pectru m 
to arbitrary ranges of k. This will be used in §1VD to 
test the effect of finite resolution and finite box size on 
the SZ power spectrum. 



IV. RESULTS 

A. Projected Maps 

Figure [l] shows a map of the density- weighted temper- 
ature for the ACDM model projected through one box 
at z = 0. Clusters of galaxies are clearly apparent as 
regions with k B T > 3 keV. The gas in filaments and 
groups can be seen to stretch between clusters and has 
temperatures in the range 0.1 < k B T < 3 keV. While 
these regions have smaller temperatures, they have a rel- 
atively large covering factor and can thus contribute con- 
siderably to the y-parameter and to the SZ fluctuations. 
This can be seen more clearly in Figure |^, which shows 
the corresponding map of the comptonization parameter. 
Clusters produce y-parameters greater than 10 -5 , while 
groups and filaments produce y-parameters in the range 
10~ 7 - 10~ 5 . Note that of the total SZ effect on the sky 
would include contributions for a number of simulation 
boxes along the line-of-sight. In such a map, the filamen- 
tary structure is less apparent as filaments are averaged 
out by projection |19| , p3| ]. A quantitative analysis of the 
contribution of groups and filaments to the SZ effect is 
presented in the following sections. 



B. Mean Comptonization Parameter 

The evolution of the density-weighted temperature for 
each of the simulations is shown on figure y. The tem- 
peratures at present are listed in table || and are quite 
similar. This is expected since all models were chosen to 
have similar as normalizations. The evolution is steeper 
for the SCDM, flatter for the OCDM model, and inter- 
mediate for the ACDM model. This is consistent with 
the different rate of growth of structure for each model. 

Also plotted on this figure is the density- weighted tem- 
perature derived from the PS formalism (Eq. [0). The 
agreement for z < 4 is good, both for the relative am- 
plitudes and for the shapes of the temperature evolution. 
At z > 4 the non-linear mass scale is not sufficiently large 
compared to the mass resolution of the simulation, so the 
temperatures are not meaningful in that regime. This is 
however not a serious limitation, as these redshifts do not 
contribute significantly to either the mean comptoniza- 
tion or the SZ fluctuations. The PS temperatures exceed 
the simulations at low redshift for all cosmological mod- 
els, since massive (high temperature) clusters, which may 
be missed in the simulations due to the effect of finite box 
size, dominate there. The PS temperatures at z = are 
listed in Table 0. 

The parameters of our ACDM model were chosen to 
coincide with that for the simulation of Cen & Ostriker 
. While the slope of our density- weighted tempera- 
ture agrees approximately with theirs for z % 3, the am- 
plitude is significantly different. They find a final temper- 
ature of about 0.9 keV, which is a factor of about 5 larger 
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than ours. This discrepancy could be due to the fact 
that their simulation include feedback from star forma- 
tion, while ours only comprise gravitational forces. It is 
however surprising that standard feedback could produce 
such a large difference. One can estimate the gravita- 
tional binding energy of virialized matter from the cosmic 
energy equation fl42| , ff3| and finds a thermal component 
of fluids to be of order 1/4 keV, consistent with our simu- 
lations. We should note, however, that the high thermal 
temperatures from feedback may be required for consis- 
tency with the X-ray background constraints p2fl . The 
reason for this discrepancy is still unknown at present, 
but should be kept in mind for the interpretation of our 
results. 

The mean comptonization parameter for each simu- 
lation was derived using equation (^|) and is listed in 
table II In all cases, y is well below the upper limit 
y < 1.5 xl0~ 5 (95% CL) set by the COBE/FIRAS instru- 
ment p8[ . The differential and cumulative redshift de- 
pendence of y are shown on figure |[ For the three mod- 
els, most of the mean SZ effect is produced at z < 2. The 
contribution from high redshift is largest for the OCDM 
and smallest for the SCDM model, again in agreement 
with the relative growth of fluctuations in each model. 

The differential and cumulative redshift dependence of 
y derived from the PS formalism are shown on this figure 
as the thin lines. The values of y from PS are also listed 
m Table || They are higher than that for the simula- 
tions by about 25% for the ACDM and OCDM models, 
are in close agreement for the SCDM model. The shapes 
of the differential curves approximately agree, although 
the PS formalism predicts more contributions from lower 
redshifts. This is can be traced to the slightly steeper 
evolution of the PS temperatures in figure |3[ and is due 
to massive nearby clusters. 



C. Power Spectrum 

As noted in jj jljj , the SZ power spectrum can be de- 
rived from the history of the temperature T p and of the 
pressure power spectrum P p (k) (Eq. 0). The evolution 
of the pressure power spectrum is shown on figure ^, for 
the ACDM simulation. The amplitude of P p (k) increases 
with redshift, while keeping an approximately similar 
shape. Perhaps more instructive is the evolution of the 



pressure bias b p (k,z) (Eq |13|), which is shown on fig- 



ure 0. For z < 1, bp is approximately independent of 
scale, in the fc-range of confidence (0.2 < k < 2 h Mpc -1 ; 
see figure |5|). The value of b p at z — and k — 0.5 h 
Mpc" 1 is listed in table fll[ For z > I, b p remains approx- 
imately constant on large scales, but is larger on small 
scales. Indeed, at early times, only a small number of 
small regions have collapsed and are thus sufficiently hot 
to contribute to the pressure. As a result, the pressure 
at high-z is more strongly biased on small scales. 

The SZ angular power spectrum derived from integrat- 



ing the pressure power spectrum along the line of sight 
(Eq. 0) is shown in Figure || for each simulation. For 
comparison, the spectrum of primary CMB anisotropics 
was computed using CMBFAST [Q , and was also plot- 
ted on this figure as the solid line. The SZ power spec- 
trum can be seen to be two order of magnitude below the 
primordial power spectrum below I < 2000, but compara- 
ble to it beyond that. Because of finite resolution and box 
size, the SZ power spectra should be interpreted as lower 
limits outside of t he /-ra nge of confidence highlighted by 
thicker lines (see §IVD). The SCDM spectrum is lower 
than that for the ACDM and OCDM. This is a conse- 
quence of the lower value of as for this model. Indeed, 
KK99 have shown that the SZ power spectrum scales as 
Ci oc f^Og/i, and is thus very sensitive on this normaliza- 
tion. This scaling relation also allows us to compare our 
results to that of the SCDM calculation of Persi et al. 
pfH . The amplitude of their power spectrum, rescaled to 
the same value of as, is within 20% of ours at I = 1000, 
while its shape is similar to ours. 

Figure || presents a comparison of the SZ power spec- 
tra derived from each of the three methods described 
in §[II. For both the simulations and the PS formal- 
ism, the SZ power spectra peak around / ~ 2000 for the 
SCDM and ACDM models, and around I ~ 5000 for the 
OCDM model. On the other hand, the constant bias 
models, which do not have a mass or k cutoff, peak at 
I ~ 10000-30000. This is explained by the fact that this 
model does not have a mass or k c utoff and has therefore 
more power on small scales. In §[VD, we will use this 



comparison to study the effect of finite resolution and 
box size of the simulations. 

For 200 < I < 2000, the simulation and PS predictions 
approximately agree for the SCDM and ACDM models. 
On the other hand, for the OCDM model, the PS pre- 
diction is a factor of 3 higher than that from the simula- 
tions in this range. This can be traced to the fact that the 
ACDM simulation yields a larger pressure bias b p (Eq. |l3] ) 
at low redshifts than the OCDM simulation. By inspect- 
ing the figure corresponding to Figure | for the OCDM 
model, we indeed noticed that more power was missing 
on small scales in this simulation. This is probably due 
to the fact that the OCDM simulation was started at a 
higher redshift (z = 100) than the other two simulations 
(z = 30). Due to truncation errors in the Laplacian and 
gradient calculations, modes with frequencies close to the 
Nyquist frequency are known to grow much more slowly 
even in the linear regime. This effect is reduced if the 
simulation is started later. 

The redshift dependence of the SZ power spectrum is 
shown in figure for the ACDM model. Most of the 
SZ fluctuations are produced at low redshifts: at I = 500, 
about 50% of the power spectrum is produced at z < 0.1, 
and about 90% at z < 0.5. The contribution of the warm 
gas in groups and filaments can be studied by examining 
figure [L2[ This figure shows the ACDM power spectrum 
measured after removing hot regions from the simulation 
volume, for several cutoff temperatures. Approximately 
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50% of the SZ power spectru m at I 
by gas with k B T < 5 keV. In §|VE 



= 500 is produced 
we show that these 
combined facts give good prospects for the removal (and 
the detection) of SZ fluctuations from CMB maps. 

The behavior of the power spectrum for / < 1000, in 
figures H and 0, agrees with the results of KK99 who 
studied the Poisson and clustering contributions sepa- 
rately. At low Vs , the SZ power spectrum is produced 
primarily by bright (low-redshift or high-temperature) 
objects, i.e., by massive clusters, and is thus dominated 
by the Poisson term. However, after subtracting bright 
clusters from the SZ map, the correlation term dominates 
the Poisson term at high redshift. Therefore, the SZ spec- 
trum on large angular, measured after subtracting bright 
spots, should trace clustering at high redshift. This in- 
teresting effect will be discussed in details elsewhere. 



D. Limitations of the Simulations 

It is important to assess the effect of the limitations 
of the simulations on these results. First, the finite reso- 
lution may lower the temperature T p , since it prevents 
s mall scale structures from collapsing. As we saw in 
§111 B , the resolution limits of the simulations correspond 
to halo masses of about 4.5x 10 12 fi m Mq. According 
to the PS formalism, the contribution to T p from halos 
with masses smaller than this limit is about 0.01 — 0.02 
keV, for z < 4 in the ACDM model. The SZ power spec- 
trum at I < 2000 is produced mainly at low redshifts, and 
is therefore little affected by this effect. Note however, 
that y, which is sensitive to small halos at high redshifts, 
is more affected. Indeed, the contribution to y by these 
halos is about 0.70 x 10 ~ 6 , assuming a gas mass fraction 
of Sl b /Q m . 

The finite box size and resolution also suppress power 
in the pressure power spectrum. As we saw in § |lll A 
and Figure [|, the simulations lack power for k < 0.1 
and k > 2 h Mpc -1 . To test the impact of this sup- 
pression, we consider the constant bias model described 



in §111 C. The total SZ power spectrum for this model 
is shown in figure |l0| as the solid line, for the ACDM 
case. This figure also shows the results of performing 
the same calculation, but after suppressing power in sev- 
eral ranges of k values. The finite box size (keeping only 
modes with k > k min = 0.1 h Mpc -1 ) reduces C; slightly 
for I < 200 and I > 20000, and thus does not have a 
very large effect. On the other hand, the finite resolution 
(k < k max = 2,5,10 h Mpc^ 1 ) reduces C/ considerably 
for I > 2000. 

The above results can be interpreted as follows. At 
a given I, the limited fc-range corresponds to a limited 
z-range, l/k max < r(z) < l/k m - m . Let us take k mm = 0.1 
and fcmax = 2 h Mpc , as relevant for the simulations. 
Then, I = 100, I = 1000 and I = 10000 correspond to 
0.02 < z < 0.4, 0.2 < z < oo and 5 < z < oo, re- 
spectively. Since most of contributions to C; come from 



z < 0.5, the finite fc m i n decreases C; only at low I, while 
&max does so over the entire I range. 

We conclude that the limitations of the simulations 
precludes us from predicting the SZ power spectrum out- 
side of the 200 < I < 2000 range. These limits can only 
be improved by using larger simulations. It is however 
worth noting that there could be more SZ power around 
/ ~ 10000. This might then be detectable by future 
interferometric CMB measurements that have angular 
resolutions around 1', intermediate in scale between the 
satellite missions and the planned millimeter experiments 
(ALMA, LMSA). 



E. Prospects for CMB Experiments 



The impact of secondary anisotropies on the upcoming 
MAP mission |4j| were studied by Refregier et al. [p5[ . 
They showed that discrete sources, gravitational lensing 
and the SZ effect were the dominant extragalactic fore- 
grounds for MAP. The dotted line on Figure || shows the 
expected noise for measuring the primary CMB power 
spectrum with the 94 GHz MAP channel, with a band av- 
erage of Al = 10. For all model considered, the SZ power 
spectrum is well below the noise. The rms y-parameter 
for the MAP 94 GHz beam (13' FWHM) is listed in ta- 
ble H for each model, from both the simulations and the 
PS formalism. The resulting rms RJ temperature fluc- 
tuations are of the order of a few fiK, compared to a 
nominal antenna noise of about 35 /iK. The SZ effect 
will therefore not be a major limitation for estimating 
cosmological parameters with MAP. 

For comparison, the residual spectrum from unde- 
tected point sources (5(94 GHz) < 2 Jy) expected us- 
ing the model of Toffolatti et al. is shown in figure || 
for the 94 GHz channel. Point sources dominate over 
the SZ effect at I > 300, but are comparable below that. 
Moreover, we have shown in figure |ll| and [l^, that about 
50% of the SZ power spectrum at / ~ 500 is produced 
at low redshifts (z < 0.1) and by clusters of galaxies 
{ksT > 5 keV). This confirms the results of Refregier et 
al., who predicted that most of the SZ effect could be re- 
moved by cross-correlating the CMB maps with existing 
X-ray cluster catalogs (eg. XBACS @, BCS @). 

Because of its limited spectral coverage, the MAP mis- 
sion will not permit a separation between the SZ effect 
and primordial anisotropies. Apart from a handful of 
clusters which will appear as point sources, it will there- 
fore be difficult to detect the SZ fluctuations directly with 
MAP. On the other hand, the future Planck surveyor mis- 
sion ^{J will cover both the positive and the negative side 
of the SZ frequency spectrum, and will thus allow a clear 
separation of the different foreground and background 
components p9|-|5"2"[ . Aghanim et al. [Q have established 
that, using such a separation, the SZ profiles of individual 
clusters can be measured down to y ~ 3 x 10- 7 . More- 
over, Hobson et al. pM have estimated that the SZ power 



G 



spectrum could be measured for 50 < I < 1000, with a 
precision per multipole of about 70%. Planck surveyor 
will therefore provide a precise measurement of the total 
SZ power spectrum. This would provide a direct, inde- 
pendent measurement of fib and of og Jl7[ , and would 
thus help breaking the degeneracies in the cosmological 
parameters estimated from primordial anisotropics alone. 
Note that this measurement might be also feasible, albeit 
with less precision, with upcoming balloon experiments 
which also have broad spectral coverage. 



F. The Missing Baryon Problem and Feedback 

The measured abundance of deuterium in low metal- 
licity systems, together with Big Bang Nucleosynthesis, 
predicts about twice as many baryons than what is ob- 
served in galaxies, stars, clusters and neutral gas pl| , p3[ . 
These "missing baryons" are likely to be in the form of 
the warm gas in groups in filaments fl42[ . This compo- 
nent is indeed difficult to observe directly since it is too 
cold to be seen in the X-ray band, and too hot to produce 
any absorption lines in the quasar spectra J54| ■ 

The SZ effect on large scale could however provides a 
unique probe of this warm gas. One can indeed imag- 
ine subtracting the detected clusters from SZ maps, and 
measuring the power spectrum of the residual SZ fluc- 
tuations, which are mainly produced by groups and fil- 
aments. For instance, if all clusters with fc^T < 3 keV 
were removed from the SZ map, the SZ power spectrum 
would drop by a factor of about 2 for I < 2000 (see fig- 
ure ]l2| ) . For the Planck Surveyor sensitivity quoted in 



§ IVE , this yields a signal-to-noise ration per multipole 
of about 1. The amplitude, if not the shape, of the resid- 
ual SZ spectrum will thus be easily detected by Planck, 
thus yielding constraints on the temperature and density 
of the missing baryons. 

In our simulations, we have only included gravitational 
forces. However, feedback from star and AGN formation 
can also significantly heat the IGM and thus affect the ob- 
served SZ effect. Valegeas & Silk S] (see also reference 
therein), have studied the energy injection produced by 
photo-ionization, supernovae, and AGN. In their model, 
AGN are the most efficient, and can heat the IGM by 
as much as 10 6 K by a redshift of a few. This results in 
a mean y-parameter of about 10 -6 , which is comparable 
to our value derived from gravitational instability alone. 
Preheating by feedback can thus increase the amplitude 
of the SZ effect by a factor of a few, and thus be eas- 
ily detected by Planck. Feedback can thus be directly 
measured as an excess in the y-parameter or in the SZ 
power spectrum, over the prediction from gravitational 
instability alone. Energy injection has a large effect on 
the gas in groups and filaments, comparatively to that in 
clusters. We may thus also detect the effects of feedback 
through the relationship between the X-ray temperature 
of groups (or their galaxy velocity dispersion) and their 



SZ temperature. These measurements would then con- 
strain the physics of energy injection. 



V. CONCLUSIONS 



We have studied the SZ effect using MMH simulations. 
Our results for the mean comptonization parameter is 
consistent with earlier work using the Press-Schechter 
formalism and hydrodynamical simulations. It is found 
to be lower than the current observational limit by about 
one order of magnitude, for all considered cosmologies. 
The SZ power spectrum is found to be comparable to 
the primary CMB power spectrum at I ~ 2000. For the 
SCDM model, our SZ power spectrum is approximately 
consistent with that derived by Persi et al. p(J, after 
rescaling for the differing values of a%. We found that 
groups and filaments (fegT < 5 keV) contribute about 
50% of the SZ power spectrum at I = 500. On these 
scales, about 50% of the SZ power spectrum is produced 
at z < 0.1 and can thus be removed using X-ray cluster 
catalogs. The SZ fluctuations are well below the instru- 
mental noise expected for the upcoming MAP mission, 
and should therefore not be a limiting factor. The SZ 
power spectrum should however be accurately measured 
by the future Planck mission. Such a measurement will 
yield an independent measurement of fib and erg , and thus 
complement the measurements of primary anisotropies 
& 

We have compared our simulation results with predic- 
tions from the PS formalism. The results from the two 
methods agree approximately, but differ in the details. 
The discrepancy could be due to the finite resolution of 
the simulations, which limit the validity of our predic- 
tions outside the 200 < I < 2000 range. We also find 
discrepancies with other numerical simulations. These 
issues can only be settled with larger simulations, and by 
a detailed comparison of different hydrodynamical codes. 
Such an effort is required for our theoretical predictions 
to match the precision with which the SZ power spectrum 
will be measured in the future. 

A promising approach to measure the SZ effect on large 
scales is to cross-correlate CMB maps with galaxy cata- 
logs ||. Most of the SZ fluctuations on MAP's angular 
scales (I < 1000) are produced at low redshifts and are 
thus correlated with tracers of the local large scale struc- 
ture. Preliminary estimates indicate that such a cross- 
correlation between MAP and the existing APM galaxy 
catalog would yield a significant detection. Of course, 
even larger signals are expected for the Planck Surveyor 
mission. This would again provide a probe of the gas dis- 
tributed not only in clusters, but also in the surrounding 
large-scale structure and therefore help solve the miss- 
ing baryon problem. Moreover, energy injection from 
star and AGN formation can produce an SZ amplitude 
in excess of our predictions, which only involve gravita- 
tional forces. The measurement of SZ fluctuations or of a 
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cross-correlation signal thus provides a measure of feed- 
back and can thus shed light on the process of galaxy 
formation. 
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FIG. 1. Density-weighted temperature for the ACDM simulation at z = 0. The temperature map was derived by projecting 
through one face of the box. 
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FIG. 2. Comptonization-parameter map for the ACDM simulation at z = 0. 
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FIG. 3. Temperature history of the gas. For each model, the density weighted temperature T p is shown for the simulations 
and for the Press-Schechter prediction. 
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FIG. 4. Differential and cumulative y-parameter versus redshift, for each simulation. The cummulative y-parameter predicted 
from the Press- Schechter formalism is shown as the thin lines. 
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FIG. 5. Power Spectrum of the DM density for the ACDM simulation. The power spectrum from the simulation (solid 
lines) is compared to that from the Peacock & Dodds (1996) fitting formula (dashed lines). The spectra, from top to bottom, 
correspond to z—Q, 0.5, 1, 2, 4, 8, and 16. 
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FIG. 6. Power spectrum of the pressure fluctuations for the ACDM model. 
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FIG. 7. Bias of the pressure b p (k, z) for the ACDM model. 
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FIG. 8. Power spectrum of the SZ effect for each model in the RJ regime, as derived from the simulations. The approximate 
range of confidence (200 < I < 2000) is highlighted by thicker lines. The power spectra outside of this range should be taken 
as lower limits. For comparison, the primary CMB power spectrum is shown for the ACDM model. The la uncertainty for the 
94GHz map channel is shown, for a band average of Al — 10. The power spectrum for the residual discrete sources (> 2Jy) for 
the 94 GHz MAP channel is also shown. 
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FIG. 9. Angular power spectrum of the SZ effect from each method: simulations, the PS formalism, and the constant bias 
model (with b p — 6). These results are shown for the RJ regime, for each cosmological models. 
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FIG. 10. Effect of resolution and finite box size on the SZ power spectrum. The ACDM power spectrum is shown for the 
constant bias model (with b p = 8.5), for several restricted ranges of k values for the pressure power spectrum P P (k, z). While the 
finite box size (0.1 < k < oo h Mpc -1 ) does not have much effect, the finite resolution of the simulations (0.1 < k < 2, 5, 10 h 
Mpc -1 ) reduce the full power spectrum (0 < k < oo) considerably outside of the approximate range 200 < I < 2000. 
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FIG. 11. Dependence of SZ power spectrum on redshift. For the ACDM model, the contribution to the SZ power spectrum 
up to a given redshift is shown. 
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FIG. 12. SZ power spectrum as a function of temperature for the ACDM model. The SZ power spectrum was calculated 
after removing regions with temperatures above the specified cutoff. 
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